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We consider in detail how the quantum-mechanical tunneling phenomenon 
occurs in a well-behaved octic potential. Our main tool will be the euclidean 
propagator just evaluated between two minima of the potential at issue. For 
such a purpose we resort to the standard semiclassical approximation which 
takes into account the fluctuations over the instantons, i.e. the finite-action 
solutions of the euclidean equation of motion. As regards the one-instanton 
approach, the functional determinant associated with the so-called stability 
equation is analyzed in terms of the asymptotic behaviour of the zero-mode. 
The conventional ratio of determinants takes as reference the harmonic oscil- 
lator whose frequency is the average of the two different frequencies derived 
from the minima of the potential involved in the computation. The second 
instanton of the model is studied in a similar way. The physical effects of 
the multi-instanton configurations are included in this context by means of 
the alternate dilute- gas approximation where the two instantons participate 
to provide us with the final expression of the propagator. 
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I. INTRODUCTION. 



The tunneling phenomenon through classically forbidden regions represents one of the 
most outstanding effects in quantum theory. Starting from the pioneering work of Polyakov 
on the subject |[|, the semiclassical treatment of the tunneling is usually presented via the 
euclidean version of the path-integral formalism. The basis of this approach relies on the 
so-called instanton calculus. In this scheme the instantons correspond to localised finite- 
action solutions of the euclidean equation of motion where the time variable is essentially 
imaginary. To sum up, one finds the appropiate classical configuration and subsequently 
evaluate the term associated with the quadratic fluctuations. In doing so the functional 
integration itself is solved by means of the gaussian scheme except for the zero-modes which 
appear in connection with the translational invariances of the system. Next one introduces 
a set of collective coordinates so that ultimately the gaussian integration is performed along 
the directions orthogonal to the zero-modes. In principle a functional determinant includes 
an infinite product of eigenvalues so that a highly divergent result is expected. However one 
can regularize these fluctuation factors by means of the ratio of determinants. 

Let us describe in brief the instanton calculus for the one- dimensional particle as can be 
found for instance in [0] . We assume that our particle moves under the action of a confining 
potential V(x) which yields a pure discrete spectrum of energy eigenvalues. In addition 
the minima of the potential satisfy V(x) = 0. If the particle is located at the initial time 
ti = —T/2 at the point X{ while one finds it when tf = T/2 at the point Xf, the functional 
version of the non-relativistic quantum mechanics allows us to write the transition amplitude 
in terms of a sum over all paths joining the world points with coordinates (— T/2,Xi) and 
(T/2,Xf). Performing the change t — > —ir, known in the literature as the Wick rotation, 
the euclidean formulation of the path-integral reads 

< x f \ exp(-HT)\ Xi >= N(T) J [dx] exp{-S e [x(r)}} (1) 

where H represents the hamiltonian of the model, the factor N(T) serves to normalize the 
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amplitude conveniently while [dx] indicates the integration over all functions which fulfil the 
boundary conditions at issue. As usual the euclidean action S e corresponds to 

,T/2 

S e = / 

J-T/2 

whenever the mass of the particle is set equal to unity. In the following we take care of the 
octic potential V(x) given by 



5 £ W> 



dr 



(2) 



V(x) = y(x 2 -l)V-4) 2 



(3) 



whose appearance can be seen in fig.l. Going to the regime in which u 2 ^> 1 the energy ba- 
rriers are high enough to decompose the physical system into a sum of independent harmonic 
oscillators. In doing so the particle can execute small oscillations around each minima of the 
potential located at x — ±1 and x = ±2. As usual the second derivative of the potential just 
evaluated at these points, i.e. V"(x = ±1) = 36 u 2 and V"(x = ±2) = 144 uj 2 characterizes 
the respective frequencies of the harmonic oscillators at issue. 
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FIG. 1. Profile of the potential. 



In terms of the discrete symmetry x — > —x which the potential V(x) enjoys, one notices 
how the four minima are non-equivalent as a whole since no connection is possible between 
the two sets represented by x = ±1 and x = ±2. Next we would like to make the description 
of the tunneling phenomenon to explain how the symmetry cannot appear spontaneously 
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broken at quantum level. In such a case the expectation value of the coordinate x evaluated 
for the ground-state is zero as corresponds to the even character of the potential V{x). 



In this section we carry things further to describe the tunneling phenomenon according 
to the tools of the euclidean version of the path-integral. At first glance it seems physically 
relevant to take into account the transition between the points Xj = 1 and Xf = 2. In doing 
so we need the explicit form of the topological configuration with Xi = 1 at tj = —T/2 
while Xf = 2 when tf — T/2. It is customarily assumed that T —>■ oo since the explicit 
solution of the problem is much more complicated for finite T. Fortunatelly the difference 
is so small that can be ignored mainly because we are interested in such a limit to obtain 
information about the energy of the first levels of the model. To get the instanton x c i(t) 
which connects the points X{ = 1 and Xf — 2 with infinite euclidean time, we can resort to 
the well-grounded Bogomol'nyi condition ||. The problem is solved just by integration of 
a first-order differential equation according to the zero-energy condition for the motion of a 
single particle under the action of —V(x). To sum up 



where the parameter r c indicates the point at which the instanton makes the jump. As usual 
equivalent solutions are obtained by means of the transformations r — ► — r and x c i{t) — > 
— Xci{ r )- It may be interesting at this point to remind how the instanton procedure allows 
only the connection between adjoint minima of the potential. In this regard we notice the 
existence of a second instanton just interpolating between Xi = — 1 and Xf = 1 (more on 
this later). Going back to x c i(t), the classical euclidean action Si associated with a such 
configuration is computed according to (|2|) so that Si = 22u/15. Now the conventional 
description of the one-instanton amplitude between Xi = 1 and Xf — 2 takes over 



II. THE ONE-INSTANTON CONTRIBUTION. 





<x f = 2\ exp(-HT)\xi = 1 >= N(T) I Det 
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exp(-Si) (5) 



Det [-(d 2 /dr 2 )+V"[x cl (r)}} 



{ Det[-(d 2 /dr 2 ) + is 2 

where we have multiplied and divided by the determinant of a generic harmonic oscillator 
of frequency v. This regularization term can be interpreted as a new amplitude given by 



<x f = 0\ exp(-H ho T)\xi = >= N(T) ^Det 



d 2 2 

■JT 2+U 



1/2 

(6) 



Fortunatelly the explicit evaluation of (|6]) is made according to the procedure exposed 
in 0. In short 

< x f = 0\exp(-H ho T)\xi = >= I -J (2 sinh z/T)~ 1/2 (7) 

Going back to the determinant of the stability equation over the instanton x c1 (t), the 
existence of a zero-mode x (r) requires the introduction of a collective coordinate. From a 
physical point of view this zero eigenvalue comes by no surprise since it reflects the trans- 
lational invariance of the system as a whole. In other words, there is one direction in the 
functional space of the second variations which results incapable of changing the action. In- 
cluding the right factor of normalization the explicit form of the zero-mode x d (t) corresponds 
to the derivative of the topological configuration, i.e. 

In addition the integral over the zero-mode itself becomes equivalent to the integration 
over the center of the instanton r c . When this change of variables is incorporated our ratio 
of determinants corresponds to || 



Det [-(d 2 /dr 2 ) + V"[x cl (r)}} 



-1/2 



\ Det[-(d 2 /dr 2 ) + u 2 } j 
I Pet* [-(d 2 /dr 2 ) + V"[x cl (r)]] \- 1/2 fs[ 

\ Det[-(d 2 /dr 2 ) + u 2 ] J V 2tt c [ ) 

where Det' stands for the so-called reduced determinant once the zero-mode has been re- 
moved. To make an explicit computation of the quotient of determinants we take advantage 



of the Gelfand-Yaglom where only the knowledge of the large-r behaviour of the classical 
solution x c i(r) is necessary ||. Starting from O and P, which represent a couple of second 
order differential operators whose eigenfunctions vanish at the boundary, the quotient of de- 
terminants is given in terms of the respective zero-energy solutions / (t) and <7 (t) according 
to 



DetO _ f (T/2) 
DetP 9o(T/2) 

whenever the eigenfunctions at issue fulfil the initial conditions 



(10) 



U-T/2) = 9o (-T /2) = 0, ^(-T/2) = ^(-T/2) = 1 (11) 
The zero-mode g {j) associated with the harmonic oscillator of frequency v corresponds 

to 

9o{j)= l - sinh[z/(r + T/2)] (12) 

so that now we only need the form of the solution / (r) associated with the topological 
configuration written in (J|). From the aforementioned x q (t) zero-mode we can write a 
second solution y {j) given by 

f T ds 

Voir) = x„(r) J q (13) 

Accordingly we may summarize the asymptotic behaviour of x (t) and y {r) as follows 

!Cexp(— 12ut) if r — > oo 
(14) 
D exp(6c<jr) if r — > — oo 

( exp(12cjr)/24cuC if r — > oo 
Voir) ~ (15) 
[ — exp(— Qlut)/12luD if r — > — oo 

where the constants C and D derive from the explicit form of the derivative of (|J). Now we 

investigate the particular solution / (r) which is the one we are really interested in. Starting 

from the linear combination of x {r) and y (^) given by 
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fo(r) = Ax {t) + By (r) (16) 
the incorporation of the initial conditions at issue leads us to 

f (r) = x (-T/2)y (r) - y (-T/2)x (r) (17) 
From this expression, which is exact, we can extract the asymptotic behaviour of f (j), 

i.e. 

fo(T/2) ~ ^ exp(3u;T) if T -> oo (18) 

Now we need to consider in detail the lowest eigenvalue of the stability equation to obtain 
the right value of the quotient of determinants. From a physical point of view we can explain 
the situation as follows: the derivative of the topological solution does not quite satisfy the 
boundary conditions for the interval (— T/2, T/2). When enforcing such a behaviour, the 
eigenstate is compressed and the energy shifted slightly upwards. In doing so the zero-mode 
x (t) is substituted for the f\(r) which corresponds to 

d 2 fx(r) 
dr 2 

whenever 



+ V"[x cl (r)}f x (r) = Xf x (r) (19) 



U-T/2) = f x {T/2) = (20) 



To lowest order in perturbation theory we obtain that 

/a(t) ~ f (r) + A ^ 



(21) 



A=0 



so that ultimately we can write that 

h(r) = f (r) + A f [x (r)y (s) - y (r)x (s)} f (s) ds (22) 

J -T/2 

The asymptotic behaviour of f°(r), x (t) and y (r), together with the condition 
/a (T/2) = provides us with the lowest eigenvalue A, namely 

A = 12luD 2 exp(—6ujT) (23) 



The final expression of the quotient of determinants requires a choice for the parameter 
v so that ultimately the frequency of the harmonic oscillator of reference is the average 
of the frequencies over the non-equivalent minima located at x — 1 and x = 2 In other 
words v = 9u. It may be interesting at this point to remark the differences with the well- 
grounded double- well model where the two minima of the potential are equivalent so that the 
aforementioned average is not necessary. However in this case the Gelfand-Yaglom method 
fixes the frequency v and subsequently the ratio of determinants is finite. Going back to the 
explicit form of x c \{r) (see (£D) we find 

C=*^, D = ^L (24) 

Armed with this information we can write the one-instanton amplitude between the 
points Xi = 1 and Xf — 2 according to 

<x f = 2\ exp(-HT)\ Xi = 1 >= 
9uA 1/2 . , „ 1/2 



— (2sinh9u;T)~ 1/i! JSx K x exp(-5' 1 ) lu dr c (25) 

7T 



where K\ stands for a numerical factor given by 



Kl = ^ (26) 

Apart from the first factor, which represents the contribution of the harmonic oscillator 
of reference, we get a transition amplitude just depending on the point r c at which the 
instanton makes the jump. According to the values of the interval T the result seems 
plausible whenever 



Sx K x exp(-Sa) w T < 1 (27) 

a nonsense condition if T is large enough. However in this regime we can accommodate 
configurations constructed of instantons and anti-instantons which mimic the behaviour of 
a trajectory strictly derived from the euclidean equation of motion. In doing so we get 
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an additional bonus since the integration over the centers of the string of instantons and 
anti-instantons is performed in a systematic way. 

To finish this section we sketch the situation as far as the second instanton of the model 
concerns. For such a purpose we take into account the one-instanton amplitude between 
Xi — — 1 and Xf — 1 which is based on the topological configuration x c2 {r) 



x c2 {r) = 2 cos 



7T 1 fe 12uJT -l 
— + — arccos 



(28) 



3 3 \e 12 ^ + 1 

whose classical euclidean action corresponds to S 2 = 76cj/5. As expected the second in- 
stanton reminds the case of the double-well potential since connects equivalent minima of 
the potential. In doing so the x (t), y (j) are symmetric so that the application of the 
Gelfand-Yaglom method is straightforward. The explicit form of the ratio of determinants 
at issue should be 

j Det'[-(d 2 /dr 2 ) + V"[x c2 (r)]] }- 1/2 
\ Det[-(d 2 /dT 2 ) + 3QuJ 2 } J 

where K 2 is given by 



S 2 K 2 oj dr c (29) 



A ' 2=12 te < 30) 



III. THE ALTERNATE DILUTE-GAS APPROXIMATION. 

As all the above calculations were carried out over a single instanton, it remains to discuss 
the complete amplitude which incorporates the physical effect of a string of instantons and 
anti-instantons along the r axis. The octic potential represents in this context a more 
complicated case since we need to include the whole scheme of non-equivalent instantons. It 
is customarily assumed that these combinations of topological solutions represent no strong 
deviations of the trajectories derived from the euclidean equation of motion without any kind 
of approximation. We wish to compute the functional integral by summing over all such 
configurations, with n instantons and anti-instantons centered at points ri, ...,r n whenever 
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T 



T 

<r n <- 



(31) 




Being narrow enough the regions where the instantons (anti-instantons) make the jump, 
the action of the proposed path is almost extremal. We can carry things further and assume 
that the action of the string of instantons and anti-instantons is given by the sum of the n 
individual actions. This scheme is well-known in the literature where it appears with the 



When evaluating the transition amplitude between X{ — 1 and Xf — 2 the total number 
n of topological configurations must be odd. As a matter of fact we can split n (odd) 
into the sum of two contributions n\ (odd) and ri2 (even) which represent the different 
possibilities associated with the existence of non-equivalent instantons. Accordingly we 
have ri\ topological configurations just interpolating between x = 1 and x = 2 oy x = — 1 
and x = —2. Of course identical situation appears in connection with ri2 where now the 
initial and final points of the trip are x — ±1. In addition we need to include a combinatorial 
factor F to count the different possibilities that we have of distributing the n instantons. 
Except for the last step which corresponds to the instanton analyzed in the previous section, 
we deal with a closed path of topological configurations starting and coming back to the 
point x = 1. As regards the instantons (anti-instantons) belonging to the first type we 
observe the formation of pairs in a systematic way due to the location of the four minima 
of the potential along the real axis. In doing so we have (ni — l)/2 + n 2 holes to fill bearing 
in mind that once the (rii — l)/2 pairs of instantons and anti-instantons are distributed no 
freedom at all remains to locate the topological configurations associated with ri2- To sum 
up 



name of dilute-gas approximation || . In addition the translational degrees of freedom yield 



an integral of the form 




(32) 




(33) 
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As regards the quadratic fluctuations over the n topological configurations at issue the 
alternate dilute-gas approximation means that 

j Pet' [-(d 2 /dr 2 ) + V"[x c1 (t)]] \ ~ 1/2 

1 Pet [- (d 2 / dr 2 ) +81u 2 ] J 



Pet' [-(d 2 /dT 2 ) + V"[x cl (r)}} 
Pet [-(d 2 /dr 2 ) + 81uj 2 } 



-1/2 



Pet' [-(d 2 /dr 2 )+V"[x c2 (r)}} 
Pet [- (d 2 /dr 2 ) + 36lu 2 } 



-1/2' 



(34) 



With all this information we can discuss the complete transition amplitude we are looking 
for once for simplicity we introduce the so-called instanton density, i.e. 

di = J Si Ki exp(-Si), i = l,2 (35) 



In other words 



<Q UJ \ 1/2 

<x f = 2\ exp(-HT) \ Xi = !>=(-) (2 sinh 9uTy 1/2 



"2 



^ [di co T]" 1 [d 2 u T] 

rti,n 2 

The sum S of (|36D can be written in terms of 

n = 2r + 1, r = 0, 1, ... 



nl 



(36) 



(37) 



n 2 = 2q, q = 0,l,... 



(38) 



so that 



(r + q)\ [d 1 uoT] 



2(r-5)+l 



[d 2 U T] 



2q 



(39) 



r =o q = (r-q)\ (2q)\ (2r + l)! 

On the other hand the best way of organizing this double sum should be the following 

i2r+l 7 



r ^ (2r + l)! 



2c/ 



(40) 



Next we can handle the sum S concerning the variable q taking advantage of |7[] 
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( r + q ) x 2q = sec[arcsin(a;/2)] cos[(2r + 1) arcsin(a;/2)] (41) 

including the transformation x — > ix to obtain that 

~ cosh[(2r + l) argsinh(s/2)] 
cosh [arg sinh ( s / 2 ) ] 

where s stands for the relative instanton density given by s = d<ijd\. In terms of a new 
variable z defined as 

z = argsinh(s/2) (43) 

it is the case that a typical value of r provides us with the final expression for S, i.e. 

~ exp[(2r + l)z] 

Now it suffices to introduce this result in (|4Tf ) so that 

= sinh[t/i^Texp(2:)] 

Next we can collect these partial results to write the final expression for the complete 
euclidean transition amplitude between the points Xi = 1 and Xf = 2, namely 

<x f = 2\ exp{-HT)\x t = 1 >= 
S^V^ 2 ln ■ u n rp\-i/2 sinhfdi u; T exp(z)] 



— (2sinh9o;T)- i/z — ^ - - ( 46 ) 

7T / V4 + S 2 

Our approach provides a new way of dealing with quantum-mechanical models which 
exhibit a more complicated structure of non-equivalent classical vacua in comparison with 
the well-grounded cases of the double-well or periodic sine-Gordon potentials where the 
equivalence of all the minima of V(x) is taken for granted |J. As regards the octic potential 
the topological solutions of the system inherit the property of non-equivalence. In any 
case the quantum fluctuations can be evaluated by means of the Gelfand-Yaglom method. 
The global effect of the multi-instanton configurations is discussed via the different kinds of 
instantons that take part to provide us with the final expression of the euclidean propagator. 
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